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Abstract 

On a Riemannian manifold (M, (?) we consider the k + 1 functions Fi, Fk,G and construct the 
vector fields that conserve Fi, Fk and dissipate G with a prescribed rate. We study the geometry 
of these vector fields and prove that they are of gradient type on regular leaves corresponding 
to Fi,...,Fk. By using these constructions we show that the cubic Morrison dissipation and the 
Landau-Lifschitz equation can be formulated in a unitary form. 
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1 Introduction 

An isolated dynamical system is a conservative system in the sense that there exist certain conserved 
quantities. Such systems are often described in a Hamilton-Poisson setting for which the energy and the 
Casimir functions are conserved quantities. In real life certain dynamical parameters are not conserved 
due to the fact that the systems are not completely isolated. Another instance when we have a dissipative 
behavior of some parameters is when control terms are added. For a large class of dissipative terms 
various mathematical formulations have been found. 

Starting from the notion of Poisson bracket, a dissipative bracket was introduced in the study of 
dissipative Hamilton-Poisson systems by M. Grmela [T^], A. Kaufman TH], and P.J. Morrison |22j- In 
[22], P.J. Morrison coined the notion of metriplectic systems which arc Hamilton-Poisson systems that 
are perturbed with a dissipation of metric type. Dissipative terms and their implications for dynamics 
have also been studied in connection to various dynamical systems derived from mathematical physics, 

see i, 0, m- 

Another type of dissipation which is called double bracket dissipation, was introduced by Brockett, 
see [TU] and [TI] . The double bracket equation is defined in a Lie algebra setting and has an important 
role in the study of various completely integrable systems, see [4], [5], [6], 0- It has been shown that 
this flow is a gradient flow on the adjoint orbits, see [6]. The metric considered is the "standard" or 
"normal" metric, see [1] and [2|. A particular example of such dissipation is given by Landau-Lifschitz 
equation. 

We show in the current paper that the cubic Morrison dissipation and the Landau-Lifschitz equation 
can be formulated in an unitary form. In Section 2 we start with the k + 1 functions Fi, Fk,G and 
construct the vector fields that conserve Fi, ...,Fk and the function G oscillates along these vector fields 
with a prescribed rate. We prove that all these control vector fields are generated by a vector field that 
we will call the standard control vector field. We apply this construction to the case when we have a 
dynamical system which has fc -I- 1 conserved quantities. We construct a perturbation which dissipates 
one of the conserved quantities and conserves the remanning k of them. 

In the Euclidean case such a dissipation was constructed in [24 using exterior algebra. In Section 
3 we extend this construction to a general Riemannian manifold and moreover, we prove that this 
dissipation is minus the standard control vector field constructed in Section 2. This generalization 
allows the study of dissipative models which have as phase space a general Riemannian manifold. 
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In Section 4 we study the geometry of the standard control vector field. In analogy with the case 
of double bracket dissipation, see [5], we prove that when restricted to regular leaves of the function 
F = (Fi, ...,Fk), the standard control vector field has a gradient formulation. On every regular leaf we 
construct a certain Riemannian metric such that the restricted standard control vector field is a gradient 
vector field with respect to this metric. First we construct a symmetric contravariant 2-tensor that is 
degenerate on the phase space. We prove that this tensor is in fact nondegenerate when restricted to 
regular leaves. The Riemannian metric on the regular leaves will be the inverse of the restriction of this 
tensor. Moreover, we prove that this metric is a conformal metric with the first fundamental form of 
the regular leaves. 

In Section 5 we prove that the standard control vector field can also be written as a scaled orthogonal 
projection on the regular leaves of the gradient vector field VG. Using this formulation, we study 
the connection between the standard control vector field generated by sets of functionally dependent 
conserved quantities. 

In Section 6 we prove that the Landau-Lifschitz equation can be regarded as a perturbed system 
with the perturbation being a standard control vector field. This perturbation was formulated before as 
a double bracket dissipation, see [19], [8], [23]. We will also show that the cubic dissipation of the rigid 
body introduced by Morrison [22] is again of the form given by a standard control vector field. Both 
systems can be described as perturbed Hamilton- Poisson systems on the Lie algebra so(3). Double 
bracket dissipation is obtained by conserving the Casimir function and dissipating the Hamiltonian 
function. Morrison dissipation is obtained by conserving the Hamiltonian function and dissipating the 
Casimir function. 

Various formulas and notations that are used throughout this paper are listed in the Appendix. 



2 Construction of the dissipation 

Having the fc + 1 functions i^i, f^, G, we construct in this section a family of vector fields that 
conserves Fi, ...,Ff^ and the function G oscillates along these vector fields with a prescribed rate. 

Let {M,g) be a Riemannian manifold and Fi, Fh,G : M — R be fc + 1 smooth functions. We 
construct a vector field u € X{M) that conserves Fi,...,Fk and dissipates G. A function Fa is a 
conserved quantity if and only if 

BF BF BF B BF B 

3x' Bx" ^ Bx°- ^ Bx^ 3x°-Bx^' 

where < • , • > is the scalar product generated by the metric g and VFa ~ g"'-' is the gradient 

vector field generated by the function F^ on the Riemannian manifold (M, g). The function G oscillates 
along the solutions of the vector field u after the rule 

dG 

— {x{t,xo)) = h{x{t,xa)), 
at 

where the rate of dissipation ft, : A/ — )• M is a smooth function and x{-,xo) is the solution of the system 
i = u{x) with the initial condition xq. We can write the above statements equivalently 

< u,VFi 

<u,VFfc>=0 (^-^^ 

< u, VG h. 

For X e M, we can decompose the vector u{x) = v{x)+w{x), where v(a;) £ SplW Fi{x), WFk{x), VG(a;)] 
is the dissipative part and w(a;) e Sp'^[\7Fi{x), ...,\i'Fk{x),S/G{x)] is the conservative part. Conse- 
quently, there exists ai{x), ...,ak{x),a{x) € M such that 

k 

v(a;) = a^{x)VF,{x) + a{x)VG{x). (2.2) 
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Ii'VFi{x), ...,'VFk{x),\/G{x) are linear independent, then locally around X S M the functions ai, afe, a 
can be chosen to be smooth and they are unique with this property. 
The algebraic system (|2.ip becomes 



f < ^F^{x),^Fl{x) > +a{x) < \/G{x),\/Fi{x) >=0 

< ^F,{x),\/Fk{x) > +a{x) < VG{x),VFk{x) >= 
I ELi "i(a^) < VG(x) > +a{x) < WG{x) , \/G{x) > = h{x) 



(2.3) 



In what follows, we will introduce several notations. If /i, fr, gi, : M — )■ M are smooth functions 
on the manifold M, we define the r x s matrix 



J.) 



< V5i,V/i > 

< Vgi,Vfr > 



< Vgs,yfr > 



(2.4) 



We solve the linear system ()2.3p for the unknowns ai, afc, a. The associated matrix is E|p,J' "'^^'gj 
and the augmented matrix is given by 



^(Fi,...,Ffc,G) 
'(Fi,...,Fi-,G) 



V 



.(Fi,...,Ffc,G) 
"(Fi,...,Ffc.G) 



\ 





h ) 



The determinant of the Gram matrix generated by the vectors Vi^i(x), Vi^fc(x), VG(x) has the 
properties that detE 



m. If detEl-^^' -'^'='^'(x) ^ 0, then rankE[^^;;;;'^';'^^(x) = rankE^^' 



linear system 



"(Fi,...,Ffc,G)\ 



'p'°Jq-^{x) > and equality holds when the vectors are linear dependent, see 

p''^'q-^{x). Consequently, the 



is compatible and according to Cramer's rule we obtain the solution 



Oii^) = (-1) 
a{x) 



i+k+l 



h(x) 



dotE 



(Fi,...,Ffc,G) 



■detE 



h{x) 



(Pi, 



Act y}'''^-- -'''k<'^h^\ 

act 2^^p_^^ p^^^,^(X) 



(Fi 

(Fi,...,FO^ 



\Fi 



'-^ ,(x) 



detEl^^'-"^''^(x) 



(2.5) 



where ^ represent the missing term. 

(Fi,...,Ffc,G), 



In the case when det E|^|' ' '^J^'g,j(a;) — 0, we will discuss the compatibility of the linear system 
If rankEj^;;:::5;^)(a;) < rankEg;;;;;^-^; (x), then rankE[^;::: JgJ (x) = rankEj^;;:::;^;g(x) 



and the linear system is compatible. If rankE 



■^'•^ ,{x) and Ej^;' 



,Ffc,G), 



{Fi, 
{Fi 



,Ffc, 
.Fk.G) 



(x) = rankE 



(Fi,...,Ffc,G) 
(Fi,...,Ffc,G) 



(x), then the matrices 



Qjix) have a common principal minor. The system is compatible if and 

(Fi,...,Ffc,G)/ 



^(Fi,...,Ffc,G) 

only if h{x) = 

On the open set fl := {x e M \ dct E|^^' ' '^^'Qj (x) ^ 0} we can use the solution found in (|2.5p and 
write the vector field in (12.21) as v = — w—rrr'Vo, where Vq S X(M) is the vector field which we 



dct S 



"(Fi,...,Ffc.G) 

will call the standard control vector field and is given by 



fe 



vo 



= ^(-l)''+'=+MetE 



j=i 



^(Fi,...,Ffc) 
(Fi,...,Fi,...,Ffc,G) 



dctE[^-::;5;vG. 



(2.6) 



For any x S M, it is straightforward to see that the set {ai{x), ...,af{x),a°{x)}, where a^{x) 



^ ^ (Fi,...,Fi,...,Ffc,G)^ ^ ^ ^ 



detE 



:Fi,...,Ffc 

(Fi 



(x) is a solution of the system (|2.3p for 



h{x) = detE 



(Fi,...,Ffc,G) 
(Fi,...,Ffc,G) 



(x). Consequently, the vector field Vq is a solution for (|2.1I) in the particular 
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case when we consider h{x) = det Yy-^p^'"" J^{x) for any x e M. For the case when we have only two 
conserved quantities, F and G for the initial system the standard control vector field has the form 

vo = - < VF, VG > Vi^ + 1 1 VF| p VG. 

As a summary of the above considerations we have the following result. 

Theorem 2.1. We have the following characterization of the control vector field u € X(M) that is a 
solution of (j2.1[) . 

/p p Q\ 

(i) Let {a; e Af I det Ylf^p^'"' p^ ^^{x) ^ 0} which is an open subset of M . Any control vector field 
u e X(f2) that satisfies (|2.ip (i.e. dissipates the function G with the rate of dissipation h) is of 
the form 

where we X{fl) withw{x) e S'p-L[VFi(2;), VF^ (x), VG(a;)]. 
(a) If the function {Fi'^-...f ■ ^ ^ R can be prolonged to a continuous function q : M ^ M., 

dot Sj^^ -Ffc,G) 

then the control vector field has the form u(x) ~ q(a;)vo(a:) + w(x), where w G X(Af) wit/i 
w(x) e 5p-L[VFi(x),...,VFfe(x),VG(x)]. 

("imJ /n i/ie particular case when h = det the control vector field that is a solution of (j2.1[) 

/las the form 

u{x) = vo(a;) + w(x), 
w/iere w € X(M),, wzi/i w(a;) G S'p-L[VFi(a;), VFfe(x), VG(x)] . 

The dissipative part of the control vector field u can be written always as v = qvo, where q is a 
smooth function defined at least on the open subset ft of M. In the paper [53] is constructed an infinite 
sequence of high-order dissipative vector fields. As a consequence of the above theorem, we obtain that 
this sequence is generated by the standard control vector field vq. 

We apply the above construction to the case when we have a dynamical system 

x^X{x), (2.7) 

where X £ Suppose (|2.7p admits -Fi, Fk, G : Af R smooth fc + 1 conserved quantities. We 

search for control vector fields u G X(Af) such that the perturbed system 

X ^ X{x)+u{x) (2.8) 



conserves Fi,...,Fk and dissipates G with a given rate. These vector fields are the solutions of (|2.II) . 
The next result gives the dissipation behavior of the function G along the solutions of the perturbed 
system (|2.8p . where u = vq + w. 



Theorem 2.2. The function G increases along the solutions of the system 

X — X{x) + \'o{x) + w(.t), 
where w G X(A/) with w{x) G Sp^[W Fi{x), Fk{x),WG{x)]. 

Proof. The standard control vector field Vq is a solution of dsn]) with h{x) ^ det i;[^J' - '^^'g](a;) > 
for any x G AI. Consequently, 

^{x{t, xo)) - h{x{t, xo)) = det s;^;;:::;t;G) (^(^^ ^o)) > o, 

where x{-,xo) is the solution of the dynamical system x — vo(x) + w(x) with the initial condition 
x(0, Xo) = xq. □ 
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3 The covariant formulation of the standard control vector field 

In the Euclidean case, a dissipation was constructed in |24| that preserves k conserved quantities of 
a dynamical system and dissipates another conserved quantity. We generalize this construction to 
a general Riemannian manifold and moreover, we prove that this dissipation is minus the standard 
control vector field Vq. 

In analogy to [53] , we introduce the following one form 

cj = (-1)" * (dFi A ... A dFk A *(dG A dFi A ... A dFk)), 

where * is the Hodge star operator associated to the Riemannian metric g. The dissipation vector field 

ttg(uj), for the Euclidean case, was introduced in [24J. Next we will prove that this dissipation is precisely 

the standard control vector field, i.e. vq = — where jig : il^{M) X{M) is associated with the 

Riemannian metric g. 

In local coordinates we have 

f)F BF BC 
dFi = —^dx''\...,dFk = — ^dx''^(^G== dx'^'+K 

By direct computation we obtain 

KdfiA...AdFfcArfG)™ , ^ f^' - f^' g°^'^..g°'°+^''°+^e;, i dx'^^^ /\.../\dx'- . 

^ {n - k - ly. Bx"^ Bx"" Bx'"'+^^ ^ ii...(fc+iifc+2...<„ 

Consequently, we have the following computation 

w (-l)"+'''*(dFi A...AdFfcA*(di^iA...ArfFfcAdG)) 



(|7.n|) 






dFk BFi 


BFk BG 




{n 


-k-iy.Bx''^ 


■■■ Qx^k dx°-^ 














(- 


-l)"+%\ BF^ 


dFk BF^ 


BFk BG 




{n 


-k-iy.Bx^i 


"'Bx''!' Bx'^^ 


"' Bx°^f= Bx^i'+i 




\9- 


\iai...akak+iPk^ 


bisi 

t-2-Pny 


„bkSk J q 

y <^Si...SkPk + 2---Pnq"'-^ 


mi 


( 




dFk dFi 


BFk BG 




in 


-k-iyBx''^ 


Bx'''' Bx°-^ 






eai. 


■■■a,k0.k + lPk + 2---Pr. 




\ '^Si...SkqPk + 2---Pn"-^ 



{n-k-iyBx^^'"Bxb>'Bx''^'"Bx^''Bx^'' + ^^ '"^ ^ i-"ai...a,a,H.i"-^ 

fe+1 
1=1 

whrrr i- — 8-Ffc dFi ^Fj^^aG_ bisi bkSk xsi ...Skai i ai 

^''^^''^ - dx>'i ■■■ dx>'k dx"! ■■■ dx-^kd^P^k+ry ■■■y "ai...akak + i"'-^ ■ 

We need to analyze the one forms wi, w^, cjfc+i. We have 

.k^, = ^ bkSkSS....Skak,.,^.k,. 

Bx'>^"'Bxf"' Bx^^'"Bx''k dx^k+l'J cn...akak + ^ 



m m mdF^ m ^s....,,. dG 
^ m mjF^ dFk ,,,, .kSk^s.-sk J^dx'^k^. 

dx^^'dx^^'dx'^^'^^y'dx'^^'^''' o"(si)---o"(sfc) g^^afc+i 

dZllQl) BFi BFk u , h , r , dFi dFk , dG , , 

Tr^-^r^9—y i>sgn{a)—^...—^)- dx""*' 

a 

= det^'}^''-'^''\dG. 
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Also 

^ ^ dx''^ dx'"'^ ^ ^ ^ 

Analogously, for i = 27fc, tj, = det S^!!''-'-^"^ ^ 

As a consequence of the above computations and of the definition of the standard control vector field 
(I2.6p . we obtain the following result. 

Theorem 3.1. On the Riemannian manifold {M,g), we have the equivalent description of the standard 
control vector field 

vo = (-l)"+itlg(*(dFi A ... A dFk A *{dG A dFi A ... A dFk))). 

4 The gradient formulation of the standard control vector field 

The standard control vector field vq is tangent to every regular leaf Lc ■— F^^(c) generated by the 
regular values of the function F :— {Fi, : M — R'^. We will endow every regular leaf Lc with a 

Riemannian metric Tc such that vo is a gradient vector field on Lc, i.e. 

In order to do this we first construct a degenerate symmetric contravariant 2-tensor T on the manifold 
M that is nondegenerate when restricted to every regular leaf Lc- The Riemannian metric Tc will be 
r — T"^ 

In what follows, we will construct the tensor T. In Riemannian geometry we can write the gradient 
vector field of the function G as VG — i^cff"^, where is the contravariant 2-tensor g~^{x) = 
gP''{x)-^ (g) gl^- constructed from the metric tensor g and i is the interior product. We recall the 
following standard results in Riemannian geometry which will be used several times in this section, 

ViJ(a) = g~\dH,a) = a(ViJ), dH{VK) =< VH,VK >, 

where a G n^{M) and H,K G e°°(M). 

We have the following contravariant 2-tensor VF^ (g) WFj : n'^{M) x rt^{M) M 

VF, ® VF, (a, 13) = a(VF,)/3(VF,). 

Lemma 4.1. For i,j G {1, k} we have the equalities 

(i) idGi^F, ® VF,) =< VG, VF, > VF, 

(U) VF, « VF, = grpgS^IE.^^ ^ _|^. 

Proof, (i) The proof is a direct computation, 

idG(VF «) VFj)(a) = VF <8) VF,- (dG,a) = dCiy Fi)a{y Fj) =< VG,\'F, > VFj(a), 

for any a G rji(Af). 

(ii) By straightforward computation we have 

BF dF 

VF «) VFj (dx^da;'^) = da;P(VF)dx«(VFj-) = g-\dF^,dxP)g-\dFj,dx'J) = ff'^V'^^^- 

□ 
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We define the symmetric contravariant 2-tensor T : Q^{M) x Q^{M) — >• R by 

T ^ (-1)^+^+1 det S; j; ^VF,+ det ^^^^g-'- (4.1) 

For proving tlie symmetry of T we notice tliat ^'',^''"'^"'"'^''1 = {^[^''' "'^'""''^l)'^ and consequently, 

( Fi F' F { F F' F ) 

det E fc; _ T,) J''"' J;' The contravariance can be deduced by tlie expressions in local 

coordinates, namely the coefficients of the symmetric 2-tensor T are 



aibi akbk„pq ^ " ^ " n ( F-t Fi 

9 ■■■9 9 Cih,---r,h,'lai...ak\^l---^t 



5Fi dFk 



k 



where we have used (|7.10p and made the notation ?7ai...afc(-P'i---^fc) ■— det 



dFi dFk 
dFi dFk 

,{Fu...,Fj,---,Fk 



Remark 4.1. For the case k — 1 and using the notation Fi = F, the expression detEj 

becomes the constant function 1 and the expression det ' becomes ||VF|p, where || • || is the 
norm generated by the Riemannian metric g. Consequently, the symmetric contravariant 2-tensor T 
has the form 

T ^ -VF ®VF +\\VF\\^g-'^. (4.2) 
In local coordinates we have the expresion 

i- -K9 g g g ) Q^a Q^b- 

For the Euclidean case, the above expression of the tensor T with F being the Hamiltonian function of 
a Hamilton- Poisson system was used in [3] . 

A 

Next, we study a few properties of the symmetric contravariant 2-tensor T that we need in what 
follows. The functions Fi, : M — R generate the following distribution on M, 

Xta„(Af) - {X e X{M) I dF,{X) = 0, s = TTfc} 

and its dual distribution 

nl^^{M) - {a e n\M) I a{VFs) = 0, s = M}. 
Proposition 4.2. We have the following results: 

(i) VF, ® WFj{a,(3) = 0, Va G f7i„„(M),V^ E n^{M). 
(ii) T(a,/3) = dets[^;;:::;^^59-i(a,/3), Va G ^\,,,{M),^P e ^\M). 
(ill) T(a,a) = detl][;;;;;:::;^;;j||a||2, Va e riLJM). 
(iv) T(dF,,/3) = 0, Vs =T7fc, V/3 e Q}{M). 
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Proof, (i) By definition we have VF^ (g) VFj(a, /3) = a{\/ Fi)l3{S/ Fj). Because a G we have the 

equality a{VFi) — 0. 

(m) From the definition of T we have that for any a E nl^„{M) and for any (3 e ft^{M), 



^(-ir+^+Meti]|;;;;;;J;;;;JJvi^^Vi^(a,/3) + detE[^;;:::;g^^^ 



= detE[^;;:::;^^;g-i(a,/3) 
(Hi) From (m) we have that for any a e f2j^„(Af) 



T(a,a) =detE[^^--^';|g-i(a,a) = detI][^^'-'-^'=|||a|i2 



(iw) From Lemma (j4.ip (i), for any s € 1, fc we obtain 



A: 



E (-1)^^^^' det sj;;;::: J;;;;; j; < v^;, v^; > v^ - + det ^-^^f, 

k 



k 



< Vi^«,VF, > (^(-l)^-+J detI]p'-'^--'^\V^^,)+detI]^-^^---!!;''|VF, 



4=1 



/ < VFi, VFi > 



= < VFs,VF, > det 



< WFi,WF, > 

< \/Fu\/Fk > 



< VFk,\7Fi > \ 

< wK^F, > 



< WFk,WFk > 
VFk 



-deti][^:;:;^:!v^; 



/<Vi^i,VFi> ... <VFfc,VFi> <VF,,VFi \ 



det 



< VFi,VFfc > 
V VFi 



<VFfc,VFfc> <VF„VFfc> 



0. 



□ 



Property (iv) of the above proposition shows that the 2-tensor T is degenerate and consequently it 
is not the inverse of any covariant metric 2-tensor. Nevertheless, the standard control vector field vq 
still behaves like a gradient vector field with respect to the degenerate symmetric contravariant 2-tensor 



Theorem 4.3. On the manifold (M,g), the standard control vector field vq is given by the following 
formula, 

vo = idG(T). 

Proof. By developing the determinants det S^''!^^' ' '^'''' after the last column in the expression 
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(|2.6p . we obtain 

k k 

= 5](-i)-'+'+^(5](-ir+'=deti]JJ;;;;J;;;;5^^< VG,V^;>)V^- + deti][^ 

j=i i=i 

k 

= E (-1)^+-'+' det sj;:};: < vg, v^^ > v^- + det vG 

- irfG(T). 

In the above proof we have used the fact that = (S^^'-'-^^' -^'h'^. □ 

The symmetric contravariant 2-tensor T is degenerate and consequently it cannot be inverted. We 
will prove that it is invertible when restricted to r2(^„(iVf). Throughout the remaining of this section 
we will assume that det ' ^ on the whole manifold M. If this assumption is not true, then 

we replace M with the open subset of regular points, i.e. {x £ M \ det s|^^' ' 'p^j(a;) ^ 0}. 

Having a Riemannian manifold (M, 5), we recall the following well known notions. One can define the 
following operators \)g : X(M) ^ rji(M), \)g{X) g{X, ■) and tJg : n^{M) -> X(M), ttg(a) := g-^{a, ■). 
The nondegeneracy of the metric tensor g implies bg — jj^ ^. For the tensor T, we can define the operator 

Ht : ^\M) ^ X(Af), BtH := T(a, •). 

We wiU prove by double inclusion the set equality '^ri^tani^^)) = '^tan{M). Indeed, for a £ 17j^„(M) 
and using Proposition 14.21 (iv), we have dFs{^T{ct)) = fiT{ct){dFs) = T{a,dFs) — which implies 
that ttT(a) e Xtari(M). For the other inclusion, let Xo £ Xtan{M) and ao := rw; Fryt'g(Xo). 

( , . . . , ) 

We will prove that tJT(ao) ~ Xq. For this, we first need to show that ao introduced above is an 
element in nl„{M). We have aoiVF^) = , ^(fI...,f,) bg(Xo)(VF,) = ^ ^^) g(Xo,Vf,) = 

fr^dFs{Xo). Because Xq £ Xta„(M), we obtain that ao{\7Fs) = 0, Vs ~ which implies 

that Q!o € By direct computation we have 

ttT(«o)(/3) =T(ao,/3)™"^ detI][;;;;;:::;^;;j5-i(ao,/3) =5"'(b<,(^o),/3) =^o(/3), 
for any ^£n^{M). 

The operator jj^ is injective on r2j^„(M). Indeed, for ai,a2 G 17f^„(M), suppose that ttrCai) = 
tJT(a2)- This is equivalent with T(ai,^ = T(a2,/3), for aU /3 e 17i(M). Using Proposition IIJ (n) , 
we obtain the equality det ^) = det I][^J;;;;;^^jg-i(a2, for all /3 € rji(M). By the 

nondegeneracy of the metric tensor g we obtain that ai = a2- 

The restricted operator (It : nl^^{M) — Xtan(-/Vf) is invertible and consequently we can define the 
inverse operator by : Xta„(A/) — >■ r2j£j„(M). From the above considerations we obtain the equality 

^TiX) - ,1^ F,M ^)^ (4-3) 

for aUX e Xt„„(M). 

Lemma 4.4. For a £ Q.^(M) and Xg £ Xtan{M) we have the equality 

T(a,bT(Xo)) = a(Xo). 
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Proof. By direct computation we have 

T(a>T(^o))™^"^dets[^;;:::;gj5-^(a,bT(Xo)) =5"'(bT(^o),a) =a(^o). 



□ 



Definition 4.1. We introduce the symmetric nondegenerate covariant 2-tensor T ^ : Xta„(Af) x 
Xta„(M) ^ e°°(M) defined by 

T-\X,Y)=T{bT{X),bT{Y)). (4.4) 

On every regular leaf Lc = F~^(c) we will construct a Riemannian metric Tc using the tensor T~^. 
Let ic ■ Lc ^ M he the canonical inclusion of the regular leaf Lc into the manifold M. We have the 
following inclusion ic,{'X,{Lc)) C Xtan{M). 

Definition 4.2. On a regular leaf Lc we define the Riemannian metric Tc '■ 'X[Lc) x X(Lc) C°°{Lc) 

Tc{X^,Y') ■.^T-\ic,X^,ic,Y^). 

The next result gives the formula for the standard control vector field by using coordinates on the 
regular leaf Lc- More precisely, we prove that Vq restricted to Lc is a gradient vector field with respect 
to the Riemannian metric r^. Moreover, we prove that this metric is a conformal metric with respect 
to the first fundamental form induced by the ambient metric g on the submanifold Lc- 

Theorem 4.5. On a regular leaf Lc we have the following characterizations. 

("i) Tc = „(Fi,...,Fi) . *c5 
(Fl Fj^) 

(a) The standard control vector field vq is a vector field in Xtan (M) . Moreover, there exists a vector 
field Vq G X(Lc) such that ic*(v3) = Vo|l^, where Vqi^^ is the restriction o/vq to the submanifold 
ic{Lc). 

(Hi) vo|L^ = ic.Vr,(Go ic). 
Proof, (i) We have the following computations 

Tc{X-,Y-) - T-\^c,X-,^cY^)^T{\>T{^c,X-).\>T{^cY')) 



detE[^;;:::;^^|g-i(bT(*c.^^)>T(*c.n) 



^ Zc,Y^) = r^,:^. Y^). 



. .^{F^.,....F^)^^"'^'- ' dets(^^---^^^o7 

°e^^(Fi,...,_Ffc) ^'^^^{F^,...,F^)°^~ 

(a) By construction, the standard control vector field vq is a solution of (|2.ip . which implies that 
Vq G XtaniM). As ic : Lc ic{Lc) C M is a diffeomorphism, there exists a vector field Vq G X(ic) 
such that ic, (v§) = vo[i,^. 

(Hi) By the definition of a gradient vector field we have 

which is equivalent with the equality 

dG{icY'=) = T-i(ze. V,,(G o ic), ic,Y'), Vr= e X{Lc). 
We have the following computations 

T-l(vo|L^,*c.n = T(bT(vo|Lj,bT(»c.n) =5"'(b9(vo|Lj>T(«c.n) = Vo|L.(bT(«c.n) 
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Consequently we have the equahty 
or equivalently 

T-i(i,,(v,^(Goi,)- vg),i,.y^) = 0, vr= e x(i,). 

And by the definition of Tc we obtain 

T,(v,,(G o I,) - vg, y=) = 0, vr= e x(Le). 

By the nondegeneracy of the Riemannian metric tensor Tc we obtain Vr^ [G o i^) — Vg or equivalently 

«c.Vr,(G O ic) = Vo|L,- □ 

For the case when we have only one conserved quantity F, i.e. k — I and Fi = F, we obtain that 

1 



IVFIP oi, 



5 The projection method formulation of the standard control 
vector field 

In this section we prove that the standard control vector field vq can be written as a scaled orthogonal 
projection on the regular leaves Lc of the gradient vector field VG. Using this formulation we study 
the connection between the standard control vector field generated by sets of functionally dependent 
conserved quantities. 

Definition 5.1. Let W he a subspace of a finite dimensional inner product space {V, < •, • >). Denote 
by the orthogonal complement ofW inV. Define Pw : V ^ V by 

Pw{v) = w, 

where v = w + u with w G W and u g . 



The linear operator Pw is called the orthogonal projection of V onto W along . 

For X e M we consider the inner product space {TxM, < •, • >), where the inner product is generated 
by the Riemannian metric g. For regular points of the function F = {Fi, ...,Ffe) : i\/ — > K.'^ we consider 
the subspace TxLc — KerDF{x). Consequently, we have TxL^ = Sp[WFi{x), ...,WFk{x)]. Indeed, for 
y e KerDFix) we have < y,VF,(x) >= g,^{x)y'g-o {x)§tix) = ^fy'ff (x) = y''§i{x) = for any 
s — l,k. 

We define the following linear operator Pt^l^ ■ T^M ^ T^M 



/ <VFi{x),VFi{x)> ... <VFkix),\/Fi{x) > <v,VFi(.t)>\ 



(^)= - ~iFu...,F,), . det 



dets|;;;:;;:^(x) 



< VFi(a;), VFfc(a;) > ... < V Fk{x) ,V Fk{x) > < v,V^^fc(x) > 
V VFi(a:) ... VFkix) v / 

(5.1) 

where v e T^M. The operator Pt^l^ is the orthogonal projection of T^M onto T^Lc along T^L^. 
Indeed, we observe that Pt^l^{^ Fi{x)) — 0, for all i — l,k, as two columns in the determinant become 
equal. Consequently, Pt^l^(u) = for any u e T^Lj: . For w e T^Lc we have 



Pt^lA^) = Tp PM '^^^ 

dets[^:;::i:J(x) 



/ < VFi(x),V^^i(x) > ... <yFk{x),yFi{x)> 0\ 

< Vi^^i(a:),Vi^^fc(a:) > ... <VFk{x),VFk{x)> 
V VFiCx) ... VFfc(x) w 



11 



Also, the standard control vector field Vq, defined by the equation (I2.6p . can be formally written as 



/ < VFi(a;),VFi(a:) > ... < WFk{x),WFi{x) > < WG{x),WFi{x) > \ 



vo(x) — det 



< VFi{x),VFk{x) > 
V VFi(a;) 



< \/Fkix),\7Fkix) > < VG{x),VFkix) > 
VFkix) VG{x) 



(5.2) 



From the above considerations we can conclude the following result. 

Theorem 5.1. The standard control vector field can be written as a scaled orthogonal projection on the 
regular leaves Lc of the gradient vector field VG, i.e. for x £ Lc, we have 



deti]|^;;:::;^^;(x)PT.L.(vG(x)). 



In what follows we study the connection between the standard control vector fields generated by two 
sets of functionally dependent conservation laws. More precisely, let 

Hi=h^{Fu...,Fu) 
Hk^hkiFi,...,Fk) ' 

where h = {hi, /ifc) : R''' — M'^ is a local diffeomorphism. By a straightforward computation we have 
the following equality, see (12. 4p . 



(_Fi,...,_Ffc) _ d{Fi, ...,Fk) ...,Ffc) -^T 



(P^'-'P") 9(xi,...,x„) 
By the same type of computation we have 

^{Hi,...,Hk) _ 9 {Hi 



d{xi,...,Xn) ' 



\Ht,...,Hk) 



d{xi,...,Xn) 



^Hk) ^^ij^^d{Hi,...,Hk)^T 



d{xi, ...,x„) 



d{hi,...,hk) d{Fi,...,Fk 



d{Fi, Fk) d{hi, hk) .t 



d{Fi,...,Fk) d{xi,...,Xn) d{xi,...,Xn) d{Fi,...,Fk) 

d{hi, hk) (F^,...,F,), ^{hl, hk) T 

d{Fi,...,Fk) (^^'■■■'^^^^d{Fi,...,Fk)' ■ 



We obtain the following equality 



(5.3) 



detS 



where c G M*^ with c = {Fi{x), Fk{x)). Consequently, on a regular leaf Lc the standard control vector 
fields Vg ^' "' ''' generated by the set of conserved quantities {Hi, ...,Hk) and respectively Vq ' ' 
generated by the set of conserved quantities {Fi, ...,Fk) differ by a constant. More precisely. 



(H^,...Mk) 



d{hi,...,hk) I . 
''''d{Fi,...,Fk)^'^ 



6 Examples 

In this section we prove that the Landau-Lifschitz equation is a perturbed system which can be put in 
the form given by the equation (12. 8p with the perturbation being a standard control vector field. This 
perturbation can also be formulated as a double bracket dissipation, see [19], [8], [23] . Also, we show 
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that the cubic dissipation of the rigid body introduced by Morrison is again of the form given by a 
standard control vector field. 



The Landau-Lifschitz equation 

One of the main objectives of the micromagnetics theory is to develop a formalism in which the 
macroscopic properties of a material can be simulated including the best approximation to the funda- 
mental atomic behavior of the material. The history of micromagnetics starts with a paper of Landau 
and Lifschitz, published in 1935, on the structure of a wall between two antiparallel domains. 

The Landau-Lifschitz equation of motion for an individual spin has the form 

M = jMxB + j^{Mx{MxB)), (6.1) 



where M is the magnetization vector, B is the magnetic field, 7 is the gyromagnetic ratio and A is the 
damping constant. Due to physical reasoning, we assume that the magnetic field '^B is of potential 
type, i.e. = for a smooth function H : R^\{{0, 0, 0)} M. Also, we suppose that ^ > 0. The 
phase-space of the problem is R^\{(0,0,0)} endowed with the Lie-Poisson bracket given by the cross 
product. An equivalent form of L-L equation the (j6.ip is given by 



7||M|p 7 



M = M xVH + — — < M, VH > M ~ -VH. (6.2) 



The unperturbed system M ~ M x VH conserves the Hamiltonian function H and the Casimir 
hmction Co = \{Ml + M| + M|). We prove that the perturbation 

^ < M, VH > M - -VH 



is the standard control vector field vq with the conserved function being F ~ ^^Cq and the dissipated 
function being G = —H. 

Using Remark 14. 11 in the case of one conserved quantity, the symmetric contravariant 2-tensor T is 
given by 

T = -VF(8) Vi^-t- ||VF|pg-i, (6.3) 
where g is the Euclidean metric on M^. Using Theorem 14.31 we obtain 
vo = idG(T) VG,Vi^ > VF+ II VFlpVG 

= <yH,.^j^>.^j^^^VH = ^^<M,VH>M~^VH. (6.4) 
Y7pf|| Y7IIMII 7 7||M||2 7 ^ ^ 

We analyze the standard control vector field vq given by the equation (16.41) restricted on a regular leaf 
Lc- In Cartesian coordinates we have the following expression of the symmetric contravariant 2-tcnsor 
field 

- ' {T(M!+Mi + Mi-Mf)4-^4-- V M.M,-|-«-|-) 



1=1 ij-=i,i^j 



Because F is a constant of motion for the perturbed system (|6.2|) , the regular leaves which are given by 
the spheres Lc F^^{c) are preserved by this perturbed dynamic. In spherical coordinates {0, ip, r) we 
have Ml = r sin6'cos(/j, M2 = r sin 6* sin 1^9, M3 = rcosO. The symmetric contravariant 2-tensor field 
T{6, (fi, r) becomes 

\ ,d d Id d , 
7r^ oO oO sm 9 dip dip 

and it is a degenerate tensor field. If we choose a sphere Lc, where r — \/Jc, then becomes a 

nondegenerate symmetric contravariant 2-tensor field on Lc- We have the coordinate expression 



\^ ,d d Id d , 
' j^c^ do do sm 9 dtp dip 
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Consequently, 



-{d0®de + SIT? 9dip (8) dip) . 



- A2 

According to Theorem 14.51 {in), the standard control vector field on the sphere Lc is a gradient vector 
field and it has the expression 

"^"'^^ ~ 72 c2 80 00 j'^c'^sm^0 dip dip' 
Remark 6.1. The induced metric on the sphere Lc is given by 



and consequently, we have 



ilg = -!-c^{d0 '^d0 + sin^ 0dp » dip) 
A 



-i*g ^ -^--5-(d6' ®d0 + SVC? 0dip ® dip). 



As stated in Theorem 14.51 (i), we obtain the equality = ||vF||^oi ^cff- 

A 

We can also write the standard control vector field Vq by using the orthogonal projector defined by 
the formula (|5.ip . More precisely, 

Pr...va . ( < v^J- > < > ) . vG - ^ < v., VG > v.. 

By using Remark 14.11 and Theorem 15.11 we obtain 

vo = ,,,\ii9 <M,VH > M - -VH. 

The Hamiltonian function H decreases along the solutions of the Landau-Lifschitz equation as G is 
an increasing function along these solutions, see Theorem 12.21 

The metriplectic dissipation of the rigid body 

The motion of a rigid body can be reduced to the translation of the center of mass and the rotation 
about the center of mass. Rotation is conveniently described in a coordinate system with the origin at 
the center of mass and the axes along the principal central axes of inertia, by Euler's equations. These 
equations can be written in the following form 

xi = (-^ - j;)x2X3+ui 

X2 = (77 - J:)xiX3 + U2 
X3 = ij;- J^)XlX2 + U3 

where Xi = JiWi, X2 = 121^2, X3 = I^lu^, are the components of the angular momentum vector, and 
Ii > I2 > I3 are the principal moments of inertia, and wi,cij2,W3 are the components of the angular 
velocity and ui,U2, U3 are the components of applied torques. The system of free rotations has the well 

known Hamilton-Poisson formulation (so(3)*, n_, where n_ is the minus Lie-Poisson structure on 

2 2 2 

so(3)* and the Hamiltonian function is given by i?(xi,X2,X3) = 2(77 + 77 + if)' [SO]- In [22], 
Morrison has introduced the following class of metriplectic dissipation for the rigid body 

u = [/i*^]VC, 



14 



where C is a Casimir of n_ and 



mix) 



/ x\ x\ X1X2 X1X3 \ 

I If + If hh^ hi:, \ 

X1X2 ^1 I ^3 X2X3 

hh 'If 71 hh 

22 

1 X1X3 X2X3 x-^ X2 I 



Next, we prove that the above control vector field u is the standard control vector field generated by 
the conserved function F — H and the dissipated function G — C. As before, by using the expression 
(14.21) of the tensor T and Theorem 14.31 we have for this case 



vo = idG(T) = {-\JH®\JH +\\\JH\\H)\IC = mVC. 

The above form of the dissipation vector field was also found in 3 . For a fixed x S R'^, the matrix 
associated with the linear operator \\V Hix)]]"^ Pt^l^ is the matrix [h^^]{x). 

In what follows we compute the standard control vector field Vq on the regular leaves Lc := H^^{c). 
A system of local coordinates (9, ip, r) for the ellipsoid Lc is given by 



a;i = r-\//i sin0 cos 93, X2 = ry' l2Smd simp, X3 — ry' cosO, 

where r = \/2c- For this choice of parameters we have 

ii„TTii9 ^ ,sin^0cos^o9 sin^0sin^(i3 cos^0, 
||Vi?||2 = 2c( ^ + ^ + ——) 

II -'2 ^3 

and the induced metric 

i*5 = 2c(/i cos^ 9 cos^ if + h cos^ 9 sin^ (p + h sin^ 9)d9 <Si d9 
+2c(/2 — Ii) sin 9 cos 9 simp cos (p{d9 ® dp + dp® d9) 
+2c sin^ 9{Ii siv? p + I2 cos^ ip)dp} (g) dp. 

For the standard Casimir function Co = 

"9 (*^1 ~^ ~^ *^ 

3) and by Theorem l4.5l (iii). the standard control 
vector field on the ellipsoid Lc is a gradient vector field which has the expression 

„ . „ 1 sin^ 03 cosmos 9 „ , 1 1 , . d 

vo|i = 2csm6'cos6'(- — )— + 2c(— - — ) sm 93COS (/s— . 

I3 I2 h o9 Ii h dp 

We observe that for the axisymmetric case Ii — /2, the dissipation takes place only in the angle 
coordinate 9, 

1 Id 
voiL, = 2csin0cos0(— - -r)^- 

7 Appendix 

Let (A/, g) be a Riemannian manifold. We will recall some standard formulas that we have used through- 
out this paper. 
Ricci symbol 

1, if ii, are distinct, and they are an even permutation of {1, r}; 
— 1, if ii, are distinct, and they are an odd permutation of {1, (7.1) 
0, otherwise. 

From the definition of the Ricci symbol, we have 

= (^1) ^ii...ikqih+i---ir~l ('''•2) 
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Generalized Kroneker (5-symbol 

1, if ii, ...,ir are distinct, and {ji, jV} is an even permutation of {ii, 
Hi ~ { *i' distinct, and{j'i, ■■■,jr} is an odd permutation of {ii, v}; (7-3) 

0, otherwise. 



Using the two definitions above, we obtain the equahty 



(7.4) 



For a fixed q we have 



5f -f ^ 5f -f -\ (7.5) 
If we do the summation after the index q from 1 to n, we obtain the formula 

and more generally 



{n-v)\ 

The determinant of an r x r matrix can be written using the formulas 



and 



ail ■.■ air 

det ( I ^ ei^,,,i^aii^a2i2---ari^ (7.8) 



ail ■•• air 

""Jl-'-Jr det I ... ... ... I ^jiii^j2i2 "'^jrir^il-'-ir' ('^•^) 



a^l ... Gi'Y^' 

Using the notation (|2.4p and developing the determinant, we obtain 



, ^(fi,...,Ffc) _ ^-Fi dF^ ,^_^ o'"^^'=detl I ('7 10^ 



The formula for the Hodge star operator * : il'^{M) — > ri"^''(M) in local coordinates is given by 



M 

(n — r)!' 



* (dx'i A ... A dx'"-) = , ff ■ ■ -ff '"^^ Cji . . . 1 . . . j-„ tfa;^"+ ^ A ... A da;^", (7.11) 



where \g\ is the determinant of the symmetric matrix associated to the Riemannian metric g. 
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